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1. Introduction 

We consider QCD with n t light flavours, n c flavours with mass m c , and nj, flavours with mass 
tnb (nic.b 3> Aqcd, >n c /nib < 1; the total number of flavours is n/ = n/ + n c + nt,). At low energies 
<C m c b it is appropriate to use the effective theory without both b and c. It has the ordinary QCD 
Lagrangian (with re-defined fields and parameters) plus l/m n cb corrections (higher-dimensional 
operators). Operators in full QCD (e.g., bilinear quark currents) are also expressed as series in 
l/ m c,fe y i a operators in the effective theory. Traditionally [1, 2, 3, 4, 5], one first decouples b 
quarks, producing an intermediate effective theory; then one decouples c quarks (Fig. 1). Power 
corrections (m c /nib) n are neglected in this approach. Its advantage is the possibility to sum leading, 
next-to-leading, etc., powers of log(mb/m c ); however, this logarithm is not really large. 

RG 

v 

m c 



Figure 1: Two-step and one-step decoupling. 

Alternatively [6], one can decouple both b and c in a single step (Fig. 1). This can be done 
at some fixed order in a s , renormalization group summation of powers of \og(mi,/m c ) is not pos- 
sible. However, now coefficients of perturbative series can be calculated as exact functions of 
m c /mb- A non-trivial dependence on m c /mf, starts from three loops, and the first power correc- 
tion (m c / mb) 2 (cc s / 7i) 3 may be of the same order as recently computed [4, 5] four-loop corrections 
(a s /7i) 4 . It is possible to combine advantages of both approaches by RG summing powers of 
\og(mb/m c ) in the leading power term [6]. 

2. Decoupling for fields and parameters of the Lagrangian 

The bare fields in the two theories are related by 

At ] = {Q l ' 2 A^\ c^ = (C) 1/2 c^, qt ] = {Q l '\t f \ (2-D 

where power corrections are neglected. Similarly, the bare parameters of the Lagrangians are 
related by 

The MS renormalized fields and parameters are related by 

a, ( "V) = A(*V) = Cj/V^M^V), (2-3) 
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where 



(2.4) 



and similarly for other decoupling coefficients. 

Let's consider, for example, the gluon field. In the full QCD (ray flavours), the bare field is 
related to the field in the on-shell scheme by A f /} = Z° s( " r) 'A& f) , where Z° s( " /} = 1/(1+ U^ f) (0) ) , 
and Tl^ f \o) contains at least one heavy-quark loop. Similarly, in the effective theory (ra/ flavours) 
A («,) = z os(«/) A (nO ; where z os(«,) = + nW(0)), and IT^O) = because it has no scale. 
Gluon propagators renormalized in the on-shell scheme in both theories are equal to the free one 
near the mass shell, and hence = A^ up to power corrections. Therefore 



7 °s(n ; ) 7 (n f ) 7 os(n,) 

£? = %^ = l+nW(0), and Wll i) = ^^ y (2-5) 

Z A Z A Z A 

Other fields are treated in the same way. The bare decoupling coefficient for light-quark masses is 
given by 

£°=Z° S [1-Z 5 (0)], (2.6) 

where the light-quark self-energy is L(p) = ]/>Lv(p 2 ) + m q oLs(p 2 ). 

In order to find the decoupling for a s , one has to consider some vertex function: Acc, Aqq, 
or AAA. They are expanded in their external momenta, and only the leading non- vanishing terms 
are retained. In the low-energy theory they get no loop corrections, and are given by the tree-level 
vertices of dimension-4 operators in the Lagrangian. In full QCD they are equal to the tree-level 
vertices times 1 where loop corrections contain at least one heavy-quark loop 1 . Then 

C(«!o /} ) = (1 +rA C - c ) 2 {Zffzt = {\ + T Aqq f (Z") 2 Z? = (1 +IW (ZD 3 . (2.7) 
The gluon self-energy up to three loops has the structure 



1 a [nf) (a [nf) 

II A (0) = - {n b m^ e +n c m; 2e ) 7>^_r(e) + P h (n b m h * e + n c m; 4e ) T F -^-F(e) 
3 7T \ 7T 



2 



+ 



{Phg+PhlTpm) (n b m b * e +n c m M 6e ) +P hh T F {n 2 b m^ e + n 2 c m^ 6e ) 



+ Pbc ( — J T F n b n c (m b0 m c0 ) 3e 



(a (n,) V 
Tf M-r(e) +.... (2.8) 



'For Aqq at the zeroth order in external momenta, there is the only structure y^. For Acc at the first order 
we shall prove this a little later. The AAA vertex at the first order in its external momenta can have, in addi- 
tion to the tree-level structure, one more structure: d aiaia2 (g ll ' li2 k^ + cycle); however, the Slavnov-Taylor identity 
{T{d^Af 1 {x),d v A v (y),d x A x (z)}} = leads to r£*$ 3 *f i£ 2 *f = 0, thus excluding this second structure. 
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The contribution involving both b and c first appears at three loops (Fig. 2) and satisfies 



Pbc {x- l )= P bc (x) , P bc (1) = 2P hh , P™ a (x^O)^ P M x 
where P b f A is the contribution of the hard region (all loop momenta ~ nib). 



,3e 



(2.9) 







Figure 2: Some contributions to the gluon self-energy with both b (thick) and c (thin) loops. 
The ghost self-energy up to three loops has the structure 

(»/) \ 2 

v0_ 
71 



II c (0) = C h (n h m b * e +n c m-J e ) C A T F ( °^-r(e) 



+ 



(C hg + C u T F ni) {ntm^ +n c mj e ) +C hh T F [jv\mjf + n 1 c m i 



cO 



+ C bc [ — ) T F n b n c (mbom c0 ) 3e 
jnboj 



CaT f 



'a M V 



The contribution with both b and c (Fig. 3) has similar properties. 




Figure 3: The contribution to the ghost self-energy with both b and c loops. 



(2.10) 



The light-quark self-energy up to three loops has the structure 



a 



(»/) 



Iy(0) = V h (n h m b * e +n c mj e ) C F T F -^F{e) 



71 



+ 



{Vhg + V M T F ni) [n h m b ° e + n c m; 6e ) + V hh T F (n 2 b m^ e + n 2 c mj e ) 



+ V bc ( — ) T F n b n c (m b0 m cQ ) 3e 



C F T F 



a 



(»/) 

sO 



71 



T(e) +•••, (2.11) 
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where V/, c comes from the diagram similar to Fig. 3; £s(0) has the same structure. 

We choose to use the ghost-gluon vertex for finding £ a . When expanded in its external mo- 
menta up to linear terms, it has the structure 



i.e., it is proportional to the tree vertex. The leftmost vertex on the ghost line singles out the 
longitudinal part of the gluon propagator: 



Therefore all corrections vanish in Landau gauge. Those with a quark loop in the leftmost gluon 
line vanish in any covariant gauge: 



--t->-*-^-*-- 






Diagrams with quark triangles produce differences: 



a 







a 





where those with a massless triangle vanish; the remaining two give [t a ,t ]. Finally, the diagram 
with the three-gluon vertex 




gives [t b ,t a ] and cancels the previous one. In result, the two-loop vertex exactly vanishes. The 
three-loop vertex vanishes in Landau gauge; it contains no diagrams with both b and c loops. 
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We used FIRE [7] for reduction of two-scale three-loop Feynman integrals. There are two 
master integrals: Fig. 4 and the same integral with a numerator. They were considered in [8, 9]. 
At m b = m c they are not independent; e expansion of the only master integral has been considered 
in [10]. 




Figure 4: The master integral. 

The decoupling relation affl = Ca a M contams 

i ("/) 

C = ! _ 1 {nhmh 2e + ncm -2 E) Tf ^_ T{£) + .... 

We re-express the right-hand side via renormalized quantities 

^r(e) = 0C t!M z ^ ) ( ai n/) (M)l e^T(l + e) l i le , 
7i Tie v / 

m b0 = Zm f) (al nf) (fl)j m b (n) , m c0 = z\„ f) (af f) (jll)J m c (fl) . 

It is convenient to use \i=m b which is defined as the root of the equation mb(m b ) = m\,. We obtain 
ajfi^ expressed via a\ ■ {fh b ), m h , and m c (mj,). Inverting the series 

K TIE \ ' 

we obtain a[ n '\pL'). It is convenient to use \x' = m c (m b ), then the result is expressed via ai (m b ) 
and 

_ m c (m h ) 
m b 

The procedure for gauge-dependent decoupling coefficients (e.g., £a) is similar but a little more 
lengthy. 

3. Decoupling for bilinear quark currents 

Until now we considered the fields and parameters of the Lagrangian. Other operators in full 
QCD also can be expressed via operators in the effective theory. Here we shall discuss bilinear 
quark currents [11]. 
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3.1 Light non-singlet currents 

First we consider light-quark currents 

& f) =ti f) r n ^o f) =Z^\a^\n))P\n), (3.1) 

where r„ = y^ 1 • • • y^"l is the antisymmetrized product of n /-matrices, and % is a non-singlet 
flavour matrix. Their anomalous dimensions are known up to three loops [12], and n appears in the 
formula as a symbolic parameter. 

These currents can be expressed via operators in the nj flavour theory: 

j^\n) = ^ V,A0i"V) + ^(lK & )- (3-2) 

It is convenient first to find the bare decoupling = 1 (neglecting power corrections), 
then 

~ zi'Vi"' V))"' 

The bare decoupling coefficient is obtained by matching on-shell matrix elements: 

2-Ov) ^ +r M) = ( C 0)-i z os W ^ +r W) > (3 4) 

where Z£ s ^ = 1 and ry^ = at zero external momenta. Therefore, 



(O =zr« /; ^r) (i+r„(a^,C j )J • (3-5) 

Decoupling of a single heavy flavour for light-quark non-singlet currents up to three loops has been 
considered in [13] (n appears in the result as a symbolic parameter). 
The vertex up to three loops has the structure 



T n = T h (n h m b « e +n c m-J e ) C F T F ( ^r(e) 



+ 



(»/) 

7T 

(r hg + r hl T Fni ) (n b m b £ £ +n c mj e ) +T hh T F {n 2 h m b £ e +n 2 c m; Q 6e ) 



+ Fbc ( — ) T F n b n c {mbomco) 
\mmj 



-3e 



(a (nf) V 
C F T F M-r(e) +-.., (3.6) 



where comes from the diagram Fig. 5a (n appears in the result as a symbolic parameter), and 

r bc ( x - 1 )=r hc (x), r foc (i)=2r /! ,, r h b ? d ( x ^o) ^r hl x 3£ . (3.7) 

For the vector current (n = 1) 

d = l, CiO*',M) = l (3-8) 
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Figure 5: Vertices of light-quark currents (the very thin line means light quark), 
to all orders 2 . For the scalar current (n = 0) 

Co° = (C°r\ Co(m',aO = CV,a0. 0.9) 

The current with n = 4 is related to n = 0, and with n = 3 — to n = 1. We can define the 
renormalized currents with 't Hooft-Veltman 75 as 

= ^jfW , W V 7^ = ^^V) , (3-10) 

because the renormalized currents 74 (/a), 73 (ju) live in 4 dimensions. The currents with "anticom- 
muting 75" are [14, 15] 

(qyfxq)^ =Z P (aM) , =Zp(«,(A0) WV^), , (3-11) 

where the finite renormalization factors Zp^ are tuned to make the anomalous dimensions of these 
currents equal to those of 70, j\ (the last anomalous dimension is 0). The anticommuting 75 does 
not influence decoupling coefficients, and therefore 

These ratios do not depend on m c /mt,. 

3.2 Light singlet currents 

Now we shall consider flavour-singlet light-quark currents (t = 1). In addition to diagrams 
discussed above, their vertex functions contain diagrams where the light quark emitted in the cur- 
rent vertex returns to the same vertex. These diagrams vanish for even n; therefore, we have to 
consider n = 1 and 3. The decoupling coefficient for the vector current is exactly 1, The anomalous 
dimension of the flavour-singlet 73 has been calculated up to three loops in [15, 16]. As compared 
to the non-singlet case, there is an additional vertex diagram Fig. 5b [17, 2]; it has the structure 

(a {nf) V 

A s C F Tjm {n b m h0 6e +n c m; 6e ) -^~r(e) . (3.13) 

2 For example, if t is diagonal with TrT = 0, then / j^dS^ over all space (dS^ is the area element of the hypersur- 
face) is a linear combination of differences of the numbers of quarks and antiquaries of various flavours; these differences 
are integer numbers, and do not change when we switch from the full QCD to the effective theory. 
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3.3 Heavy currents 

Finally, we consider ft-quark currents (results for c-quark currents can be obtained by the 
obvious substitution). They can produce flavour-singlet light-quark currents. The number of y 
matrices on the light-quark line is always odd, hence we should consider the currents with n = 1 
and 3. The vector current simply becomes in the effective theory, up to power corrections. The 
current with n = 3 becomes the flavour-singlet light-quark current: 



j b o = bo7 [a r P f ] bo., 



Jbb 



(3.14) 



where £° starts from two loops (Fig. 6a). Decoupling of the single heavy flavour (b) has been 
calculated up to three loops in [16, 17]. When decoupling both b and c, 



A + \ A g +AiT F m +A h T F n b +A bc 
xC F T F n b (°^r(s)m h0 2e ) 2 , 



T F n t 



(%s0 
71 



r(e)m 



-It 
bo 



+ ••• 



where A bc comes from Fig. 6b, and 



A bc (l)=A h , A\f(x^Q)=A 



i hard 
l bc 



■ A s x 



-6e 



(3.15) 



(3.16) 



(in the hard region all loop momenta are of the order of the heaviest quark mass). Combining the re- 
sult for A bc with the previously known single-scale results, we obtain the renormalized decoupling 
coefficient. 



A 




Figure 6: Vertices of fr-quark currents. 



4. Conclusion 

We have considered simultaneous decoupling of b and c quarks with three-loop accuracy for 
the fields and parameters of QCD Lagrangian [6] (a s , light-quark masses m q , gauge parameter a, 
fields A, c, q), as well as for bilinear quark currents [11] (flavour non-singlet and singlet light- 
quark currents, heavy-quark currents). The coefficients of the leading power correction to the 
three-loop term (m c / m b ) 2 (a s / n) 3 happen to be of order 1, while the coefficients of the four-loop 
corrections [4, 5] (a s / n) 4 are numerically large. Therefore, our power corrections are numerically 
negligible (but this conclusion could not be obtained without calculating these power corrections). 

Acknowledgements. A.G. is grateful to DESY Zeuthen and TTP Karlsruhe for financial 
support which allowed to attend the conference. This work was supported by the BMBF through 
Grant No. 05H09VKE. 
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